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The unusual property of negative thermal expansion is of fundamental interest and may be 
used to fabricate composites with zero or other controlled thermal expansion values. Here 
we report that colossal negative thermal expansion (defined as linear expansion < -10~ 4 K~ 1 
over a temperature range -100 K) is accessible in perovskite oxides showing charge-transfer 
transitions. BiNi0 3 shows a 2.6% volume reduction under pressure due to a Bi/Ni charge transfer 
that is shifted to ambient pressure through lanthanum substitution for Bi. Changing proportions 
of coexisting low- and high-temperature phases leads to smooth volume shrinkage on heating. 
The crystallographic linear expansion coefficient for Bi 095 La 0 05 NiO 3 is -137*10" 6 K -1 and a value 
of -82x10" 6 K -1 is observed between 320 and 380 K from a dilatometric measurement on a 
ceramic pellet. Colossal negative thermal expansion materials operating at ambient conditions 
may also be accessible through metal-insulator transitions driven by other phenomena such as 
ferroelectric orders. 
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Useful negative thermal expansion (NTE) materials for zero 
or controlled expansion composites should show a smooth 
contraction while heating through a wide temperature 
range 13 . Volume shrinkages in framework-type materials such as 
ZrW 2 O g (ref. 4) and Cd(CN) 2 (ref. 5) can extend over very wide 
ranges, for example, ZrW 2 0 8 shows a near- continuous volume 
decrease of AV/V= -2.0% over 0.3-1050 K 6 , equivalent to an aver- 
age linear expansion of oc L = (1/L)(AL/AT) = - 7xlO~ 6 K _1 where the 
linear dimension L here is the cubic unit cell parameter. NTE can 
also result from transitions between different electronic or magnetic 
states strongly coupled to the lattice, giving large negative expan- 
sions down to a previous record dilatometric value of - 25xlO" 6 K _1 
for (Mn 096 Fe 004 ) 3 (Zn 05 Ge 05 )N at 316-386 K (ref. 7). Some rare earth 
compounds show similar volume shrinkages of a few % at first order 
valence transitions, for example, SmSe (ref. 8), Sm 2 75 C 60 (ref. 9), and 
Yb 8 Ge 3 Sb 5 (ref. 10), but the volume usually collapses over a narrow 
interval at low temperatures with strong hysteresis. 

Volume shrinkages resulting from intermetallic charge transfers 
were recently found near room temperature in two perovskites syn- 
thesized at high pressure, LaCu 3 Fe 4 0 12 (ref. 11) and BiNi0 3 (ref. 12). 
Although the electronic transitions are first order, the accompanying 
volume change has an appreciable thermal width for the latter com- 
pound suggesting that inhomogeneities or chemical substitutions 
could be used to make the transition broad and quasi-continuous, 
leading to large NTE over an appreciable temperature range. We 
investigate here the possibility of preparing colossal NTE (CNTE) 
materials, defined as having a smooth change of AL/L < - 1.0% over a 
temperature interval of AT-100K (and hence o^K - lOOxlO^K" 1 ), 
and we show that doped- BiNi0 3 ceramics may enable this concept 
be realized at ambient conditions. 

Results 

High-pressure and -temperature PND and XAS studies of 
BiNiO s . Intermetallic charge transfer in the perovskite BiNiO s was 
discovered by application of pressure. At ambient pressure, BiNi0 3 
has a triclinically distorted crystal structure with the unusual 
valence distribution Bi 3+ 0 5 Bi 5+ 0 5 Ni 2+ O 3 and shows charge ordering 
of the Bi valences 13 . This phase is insulating and ferrimagnetic, with 
a near-antiferromagnetic alignment of S = 1 Ni 2+ spins up to a Curie 
temperature of 300 K (ref. 14). Previous powder neutron diffraction 
(PND) and X-ray absorption spectroscopy (XAS) studies revealed 
a pressure-induced melting of the Bi-charge disproportionation at 
3-4 GPa and a simultaneous Ni to Bi charge transfer accompanied 
by a structural change to the orthorhombic GdFe0 3 type perovskite 
superstructure with valence distribution Bi 3+ Ni 3+ 0 3 (refs 12,15). 
This transition is accompanied by a broad insulator-to- metal 
transition, and a substantial volume decrease of 4.5% was observed 
between 2.1 and 4.7 GPa. A temperature-induced insulator- 
metal transition was reported for La- substituted Bi 0 95 La 0 05 NiO 3 at 
ambient pressure 16 . This transition was attributed to the hole doping 
to Ni 2+ ions from photoemission and XAS studies 1718 whereas the 
unit cell volume shrinkage of 2.9% found on heating suggested that 
substantial charge transfer occurs. 

We have reinvestigated the first order structural transition 
accompanying the charge transfer in BiNi0 3 using small pressure 
steps to estimate the transition width (Methods). Figure la shows 
part of the time -of- flight PND data at selected pressures up to 
5.5 GPa on pressurizing. The drastic change around 3.5 GPa shows 
the transition from the triclinic Bi 3+ 0 5 Bi 5+ 0 5 Ni 2+ O 3 structure (in the 
space group P-l) to the orthorhombic Bi 3+ Ni 3+ 0 3 high-pressure 
phase (in the space group Pbnm). Lattice parameters and unit cell 
volumes derived by fitting the data are shown in Figure lb,c. These 
show that all the lattice parameters decrease across the transi- 
tion, although the magnitudes of the changes differ, and the unit 
cell volume decreases by 2.5 %. This large change results from the 
dominant contraction of the Ni-O perovskite framework as Ni 2+ 
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Figure 1 1 Structural evolution of BiNi0 3 under pressure at room 
temperature, (a) Selected time-of-flight PND data at variable pressures 
showing the change in diffraction intensities around the 3.5 GPa insulator- 
metal transition, also a peak from the Pb pressure standard (shaded), 
(b) and (c) Pressure dependence of the lattice constants (b) and the unit 
cell volume (c) of BiNi0 3 . Insets show the low pressure and temperature 
(LPT, Bi 3+ 05 Bi 5+ 05 Ni 2+ O 3 ) and high-pressure and temperature (HPT, 
Bi 3+ Ni 3+ 0 3 ) crystal structures, with purple and yellow spheres respectively 
corresponding to Bi 3+ and Bi 5+ states in the LPT structure. Blue and 
red symbols on (b) and (c) correspond to the LPT and HPT phases 
respectively. The circles, squares and triangles in (b) respectively stand for 
a, b and c/V2 lattice constants. The open point in (c) shows the weighted 
average volume at 3.3 GPa. 

is oxidized to the smaller Ni 3+ at the transition, which outweighs 
the lattice expanding effects of reducing Bi 5+ to Bi 3+ and increases 
in the Ni-O-Ni angles. The structural transformation was found to 
occur from 3.2 to 3.7 GPa on compression so the transition width 
is taken as Ap CT ~0.5GPa, which is consistent with the weighted 
average volume of the two coexisting phases observed at 3.3 GPa as 
shown on Figure lc. 

The charge-transfer transition in BiNi0 3 has a negative dT/dp 
slope, which enables the same structural change and accompa- 
nying volume reduction to be observed by heating at a moderate 
pressure. Figure 2 shows the temperature dependence of the lattice 
parameters and unit cell volume of BiNi0 3 on heating at a pressure 
of 1.8 GPa, with the transition occurring around 500 K. Structure 
refinements from neutron data, XAS, and resistivity measurements, 
shown as Supplementary Figures SI -S3 and Supplementary Table 
SI, confirm that the high-pressure and -temperature (HPT) state 
is the same orthorhombic and metallic Bi 3+ Ni 3+ 0 3 phase as gen- 
erated by pressurization at ambient temperature, and none of the 
alternative charge states proposed elsewhere 12 were observed. The 
volume drop of 3.4% observed by heating in Figure 2b has a similar 
value to the pressure-induced change in Figure lb. Coexistence of 
the low-pressure and -temperature (LPT) Bi 3+ 05 Bi 5+ 05 Ni 2+ O 3 phase 
and the HPT Bi 3+ Ni 3+ 0 3 form is observed at 490 and 510 K, which 
enables the transition width of the volume decrease to be estimated 
as AT CT = 80 K by fitting a straight line through the two weighted 
average volumes as shown in Figure 2b. The volume (bulk) ther- 
mal expansion coefficient is oc v = (l/y)(AV/AT CT ) = -282xlO~ 6 K _1 
corresponding to (\- -94xlO" 6 K _1 over an 80 K range which 
approaches the CNTE paradigm. 
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Figure 2 | Structural evolution of BiNi0 3 at elevated temperatures under 
a constant pressure of 1.8 GPa. (a) and (b) Temperature dependence of 
the lattice constants (a) and the unit cell volume (b) of BiNi0 3 at 1.8 GPa 
showing the LPT (Bi 3+ 05 Bi 5+ 05 Ni 2+ O 3 ) to HPT (Bi 3+ Ni 3+ 0 3 ) phase transition 
around 500 K. Blue and red symbols correspond to the LPT and HPT 
phases respectively. The circles, squares and triangles in (a) respectively 
stand for a, b and c/V2 lattice constants. The weighted average volumes 
shown for the two points in the coexistence region in (b) are used to 
estimate the transition width by extrapolation from the 560 K point. 

Our neutron diffraction and XAS results have been used to con- 
struct the P-T phase diagram for BiNi0 3 shown in Figure 3. BiNi0 3 
decomposes above 500 K at ambient pressure, but is stabilized up to 
at least 565K at 1.8 GPa (and to -1,300 K at 6GPa under synthesis 
conditions). The boundary between the LPT and HPT phases has 
slope dT CT /dp= - 140KGPa _1 . The 2.5-3.4% volume contraction 
occurs over observed transition pressure and temperature widths of 
0.5 GPa and 80 K respectively in BiNi0 3 , which are consistent with 
the transition slope as - AT CT /Ap CT = - 160 KGPa _1 ~dT CT /dp. 

Powder XRD and dilatometry of Bi 0 95 La 0 05 NiO 3 at ambient pres- 
sure. The large AT CT of BiNi0 3 shows that CNTE is feasible but 
the transition is only observed above pressures of 1.5 GPa in pure 
BiNi0 3 . However, chemical substitutions for Bi may be used to sup- 
press the charge disproportionation in the Bi 3+ 0 5 Bi 5+ 0 5 Ni 2+ O 3 phase 
and thereby shift the charge-transfer transition to near ambient 
conditions. This is conveniently achieved by a partial substitution of 
La for Bi that suppresses the characteristic insulator to metal tran- 
sition accompanying charge transfer to around 350 K at ambient 
pressure in Bi 0 95 La 0 05 NiO 3 (ref. 16). The thermal expansion of this 
material has been explored using both X-ray diffraction (XRD) and 
dilatometry measurements to discover whether the La- substitution 
broadens the transition enabling wide-temperature range CNTE 
behaviour to be obtained. 

X-ray diffraction measurements on Bi 0 95 La 0 05 NiO 3 in Figure 4a 
show the 2.9% volume shrinkage, which has a similar magnitude 
to that observed in undoped BiNi0 3 under a pressure of 1.8 GPa. 
Coexistence of the low- and high -temperature phases is observed 
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Figure 3 | Pressure-temperature phase diagram of BiNi0 3 determined by 
PND and XAS studies. Circles and squares show PND and XAS data, and 
blue and red symbols correspond to the LPT and HPT phases, respectively. 
The width of the transition between LPT (Bi 3+ 05 Bi 5+ 05 Ni 2+ O 3 ) and HPT 
(Bi 3+ Ni 3+ 0 3 ) phases is indicated by the shaded area. 

at three points in the transition region, and a linear fit to the 
weighted average volumes is used to obtain the transition width of 
AT CT = 70 K, which is comparable to the 80 K estimate for BiNi0 3 
under pressure. The crystallographic volume thermal expansion 
coefficient between 300 and 370K is 0^= -413xlO" 6 K- 1 and the 
linear coefficient is a L =-137xlO" 6 K" 1 , showing that CNTE mag- 
nitudes are also observable in Bi lx La x Ni0 3 . Crystallography pre- 
dicts the upper limit of the magnitude of thermal expansion as the 
formation of pores and other microstructural defects can lessen the 
effect in bulk ceramics. Dilatometric measurements on a polycrys- 
talline ceramic of Bi 0 95 La 0 05 NiO 3 prepared at high pressure were 
made during heating and cooling cycles as shown in Figure 4b. 
The strain AL/L(400K) increases with increasing temperature up 
to 270 K indicating the normal positive thermal expansion, but 
decreases above 270 K. The average observed oc L between 270 and 
400 K is -49xlO" 6 K _1 and the maximum negative slope between 
320 and 380 K corresponds to a linear thermal expansion coef- 
ficient of -82xlO" 6 K _1 . This is more than three times larger than 
the previous maximum NTE value by dilatometry of-25xlO" 6 K _1 , 
reported for the linear expansion of an anti-perovskite manganese 
nitride (Mn 096 Fe 004 ) 3 (Zn 05 Ge 05 )N ceramic over a similar tempera- 
ture range (3 16-386 K) (ref. 7). Furthermore, the derived bulk ther- 
mal expansion is only 60% of the predicted crystallographic value, 
suggesting that Bi 0 95 La 0 05 NiO 3 ceramics showing CNTE should be 
accessible through further materials processing, for example, by 
using smaller grains to minimize pore formation. The hysteresis of 
up to 20 K between warming and cooling cycles would have to be 
further reduced for practical applications. 

Discussion 

BiNi0 3 -based ceramics have several advantages as negative expan- 
sion materials, in addition to colossal <\ values. They are durable, 
and our samples have proved air- and moisture -stable over periods 
of several years to date. Our study of Bi l x La x Ni0 3 has shown that the 
temperature and width of the charge transfer transition are easily 
tuned around room temperature, and many further substituents for 
Bi are available, such as other lanthanides. The volume contraction 
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Figure 4 | Thermal measurements for Bi 0 95 La 0 05 NiO 3 at ambient 
pressure showing thermal expansion coefficient values over the linear 
regions indicated, (a) Temperature dependence of the unit cell volume 
on heating with LPT and HPT data drawn as open and closed circles, and 
crosses showing the weighted average volume in the transition region, 
(b) The dilatometric linear thermal expansion of Bi 095 La 005 NiO 3 on heating 
and cooling showing a 20 K hysteresis. 

also coincides with an insulator to metal transition that could be used 
as an additional response in smart materials applications. Chemical 
tuning is also needed to suppress hysteresis. However, if hysteresis 
is completely suppressed, the phase transition will no longer allow 
coexistence of both phases, and the overall NTE effect will be lost. 
The possible solution is to drive the transition to second order, 
as was done for Mn 3 (Cu l x Ge x )N (ref. 7). The greatest limitation 
is that high pressures of ~6GPa are presently required to synthe- 
size BiNi0 3 -based materials, but non- equilibrium methods may be 
developed for bulk synthesis, as was done for diamond growth. 

In conclusion, this study shows that charge-transfer transitions 
in transition metal oxides may enable CNTE materials operating at 
ambient conditions to be realized. Although the Bi 3+ 0 5 Bi 5+ 0 5 Ni 2+ O 3 
to Bi 3+ Ni 3+ 0 3 transition in BiNi0 3 is first order, the volume shrink- 
age occurs over a reproducible width of -80 K that is stabilized to 
ambient conditions by partial substitution of La 3+ for Bi 3+ . The lin- 
ear thermal expansion coefficient of -82xlO~ 6 K _1 measured over 
a 60 K interval in Bi 0 95 La 0 05 NiO 3 by dilatometry is more than three 
times that of other materials showing NTE near room temperature, 
and further substitutions and processing may enable the crystallo- 
graphically observed magnitude of < - 100xlO~ 6 K _1 to be observed 
with suppression of thermal hysteresis. We note that other transitions 
involving the electronic degrees of freedom in transition metal oxides 
can also give rise to large volume shrinkages, for example, insulator to 
metal transitions induced by pressure in PbV0 3 (ref. 19) and BiCo0 3 
(ref. 20) lead to the loss of off-centre ferroelectric displacements, 
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and a spin state transition in the latter material, with > 10% 
volume reductions. Chemical substitutions to shift these transitions 
to ambient conditions and broaden them into continuous changes 
could result in other CNTE materials. 

Methods 

Sample preparation. Polycrystalline BiNi0 3 was prepared at 6GPa and 1,000 °C 
with a cubic anvil type high-pressure apparatus, as reported previously 13 . 

PND measurements. PND patterns at high-pressure, high-temperature condi- 
tions up to 7.7 GPa and 560 K were taken with a Paris -Edinburgh press at the Pearl 
beamline of the ISIS facility, UK. Pb was used as a pressure marker in ambient 
temperature experiments. The temperature variation of an 80 mg sample was 
explored at 1.8 GPa pressure using a variable-temperature insert. MgO was used as 
a pressure marker and the temperature was determined by the neutron resonances 
of Hf and Ta 21 . Data were collected for approximately 1 h at each temperature, with 
rapid heating between successive temperatures. Rietveld fits to the data were used 
to fit the lattice parameters and, in the coexistence regions, phase fractions that 
were used to calculate the weighted average of cell volumes. Representative high- 
pressure neutron diffraction results are shown in Supplementary Figure SI. 

XAS measurements. The XAS measurements on BiNiO s at the Ni-K edge under 
various pressures and temperatures were carried out using synchrotron radiation at 
the BL39XU beamline of SPring-8, Japan. A diamond anvil cell having anvils with 
a culet diameter of 0.6 mm was used to apply pressures up to 5 GPa. The powdered 
sample was mixed with NaCl as a solid pressure medium. The pressure inside the 
DAC was determined from the NaCl lattice constant. Sample heating was achieved 
with a resistance heater placed around the sample. 

Synchrotron X-ray powder diffraction measurements. The temperature 
dependence of the lattice constant of Bi 0 95 La 0 05 NiO 3 was determined from synchro- 
tron X-ray powder diffraction data collected at beamline BL02B2 of SPring-8 with 
wavelength X = 0.42084 A. Less than 1 mg of sample powder was loaded into a glass 
capillary with an inner diameter of 0.1 mm. The sample temperature was controlled 
by blowing cold or warm nitrogen gas. The heating rate was 60 K per min, and the 
temperature was stabilized for 5 min before each measurement. 

Dilatometrc measurements. The linear thermal expansion of a 3 mm diameter 
and 1 mm high pellet was measured by dilatometry with a KYOWA type KFL strain 
gauge between 100 and 400 K. The temperature dependence of the resistance of the 
strain gauge was measured with a Quantum Design PPMS at heating/cooling rates 
of 1 K per min. Cu and Si0 2 were used as reference materials. 
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